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Magnetoresistance (MR) has attracted tremendous attention for pos¬ 
sible technological applications. Understanding the role of magnetism 
in manipulating MR may in turn steer the searching for new applicable 
MR materials. Here we show that antiferromagnetic (AFM) GdSi metal 
displays an anisotropic positive MR value (PMRV), up to ~ 415%, accom¬ 
panied by a large negative thermal volume expansion (NTVE). Around 
Tii the PMRV translates to negative, down to ~ -10.5%. Their theory¬ 
breaking magnetic-field dependencies [PMRV: dominantly linear; nega¬ 
tive MR value (NMRV): quadratic] and the unusual NTVE indicate that 
PMRV is induced by the formation of magnetic polarons in 5d bands, 
whereas NMRV is possibly due to abated electron-spin scattering result¬ 
ing from magnetic-field-aligned local 4/ spins. Our results may open up a 
new avenue of searching for giant MR materials by suppressing the AFM 
transition temperature, opposite the case in manganites, and provide a 
promising approach to novel magnetic and electric devices. 
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Magnetoresistance (MR), a change in electrical resistivity when an external magnetic 
held {hqH) applied, occurs in metals, inorganic and organic semiconductors, and particularly 
close to an intermediate regime of the transformations of charge (insulator, metal) order and 
spin [paramagnetic (PM), ferromagnetic (FM)] order in thin manganites’ hlms as a colossal 
negative MR value (NMRV) [IH5]. The colossal MR (CMR) in manganites is accompanied 
by a shift of the transition temperature to a higher value by applied magnetic held so that 
a sharp peak appears in the MR values near the transition, e.g., the NMRV reaches ~ -95% 
at po-R = 15 T near Tc = 240 K in Lao.ssSro.isMnOa [6]. Although the double exchange 
interaction can qualitatively explain the CMR ehect based only on the charge and spin 
degrees of freedom, the exact picture still remains elusive. 

We investigate here the intermetallic GdSi that crystallizes in the FeB-type structure 
{Pnma) and orders antiferromagnetically below ~ 55 K [5l 0^. The vanished orbital mo¬ 
mentum {L = 0) naturally removes the crystal electric held (CEF) as well as its perturbation 
on magnetic interactions, which renders GdSi ideal for studying pure spin-magnetism-tuned 
transport. Meanwhile, there keep solely two main sources for the magnetovolume (MV) 
ehect. One contribution is from the volume dependence of the Ruderman-Kittel-Kasuya- 
Yosida (RKKY) exchanges, which normally produces a negative MV (NMV) ehect. The 
other is due to the spin-polarized itinerant moments in the conduction bands, resulting in a 
positive MV (PMV) ehect. Therefore, the volume variation (NMV or PMV) with temper¬ 
ature in a magnetic state could be useful in determining the nature of its magnetic origin. 
Even after nearly 50 years of research on GdSi [SUM, to our knowledge, there is little data 
on its temperature-dependent structural modihcations and MR property. Here we report on 
detailed temperature-dependent powder dihraction, angular-dependent magnetic character¬ 
istic and transport studies of GdSi single crystals. We hnd anomalous anisotropic giant MR 
and spontaneous magnetostriction (MS) ehects, in particular, an antiferromagnetic (AFM)- 
driven negative thermal volume expansion (NTVE) and a nontrivial positive to negative MR 
transition, which can be well understood by combining the magnetic tunnel of conduction 
electrons and the concept of magnetic polarons Ha, i.e., local short-range FM spin regimes 

na. 


Results 

Structural studies. The x-ray powder dihraction analysis (Fig. 1) clearly shows three 
distinct structural regimes (I, II and III). Upon cooling, the rehned (Re) a, b, c and V shrink 
almost linearly before Ti, followed by a slower decrease until sharp turns at Tm, an onset 
temperature of the AFM transition. Below Tjn, the increases of a and b and the decrease 
of c ultimately result in an unusual NTVE (i.e., PMV) in the unit-cell volume F, which 
is quite useful in fabricating applicable materials with controlled thermal expansion values 
PSIIII]. Obviously broadening in the nuclear Bragg peaks is present in the I regime (see 
Supplementary Fig. Sib), which is attributed to the magnetoelastic ehect (the coupling 
between magnetic moments and lattice strains). The strain distribution patterns are ex¬ 
tracted and shown in Supplementary Fig. S2. The variations in a, b, c and V below Tm 
imply that magnetic anisotropy (MA) (which is consistent with Supplementary Fig. S3a and 
uncommon for the F-state Gd-compounds), spontaneous PMV and anisotropic MS ehects 
exist in GdSi. 

The GEF is mainly responsible for the giant MS ehect in rare-earth (RE) compounds. 
This ehect is thus expected to be negligible in GdSi. However, below Ti, structural param¬ 
eters shown in Fig. la obviously deviate from the the oretical estimates (solid lines) by the 


Griineisen (Gr) law, e.g.. 


= 4.94(6) X 10 
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FIG. 1. Temperature-dependent structural parameters, (a) Anisotropic character of the 
lattice-constants, a, b and c, variation, (b) Anomalous unit-cell volume, V, expansion with tem¬ 
perature in the Pnma symmetry. The solid lines are theoretical estimates of the temperature- 
dependent structural parameters using the Griineisen model with Debye temperature of 0 d = 340 
K, which is the same as reported in Ref. [5j. Upon warming, two appreciable anomalies display in 
the structural parameters at respective temperatures of Tjn ss 54.3 K (0.06 T) and Ti « 100 K. 
Error bars in (a) and (b) are the standard deviation obtained from the Fullprof refinements. 


-2.60(8) X 10“^, denoting large anisotropic spontaneous MS effects. The formation of long- 
range-ordered (LRO) AFM state is a process of the growth of sublattice FM domains. The 
enlargement of FM domain volumes with decreasing temperature may accumulate strains 
on the domain walls, which is the microscopic mechanism for the magnetic-held-induced 
MS effect in ferromagnets. Therefore, including the effect of the molecular held of one Gd- 
AFM-sublattice on the other is indispensable to understand the spontaneous MS effect in 
GdSi. According to the Stoner model for itinerant magnetic electrons, the positive magnetic 
pressure Pm associated with the magnetic ordering in a band is proportional to 
where D, V and M represent the electronic density of states at the Fermi energy, the vol¬ 
ume and the magnetic moment, respectively. The spontaneous PMV effect (i.e., NTVE), 
e.g., = 1.04(2) X 10“'*, is therefore attributed mainly to the increases of D (cor¬ 

responding to the pronounced decrease of p below Tu in Fig. 2a) and the induced itinerate 
spin-moments in conduction bands. Similar MS and MV effects were also reported in GdsNi 
[18], where, however, they are ascribed to the itinerant character of the Ni 3d bands. 

Resistivity measurements. As depicted in Fig. 2a, the zero-held-cooling (ZFG) elec¬ 
trical resistivity at 0 T with current / along the a and c axes decreases linearly due to the 
weakened thermal excitations upon cooling until around Tu (vertical bar), where a hump 
with negative slop appears along only the c axis, probably attributed to the magnetic su¬ 
perzone effect [19] as a consequence of the AFM ordering. Below Tn, they decrease steeply 
like a reasonable AFM metal. The resistivity between 10 K and 40 K [6] can be well £t to 
p{T) = po + kT^, shown as dashed lines. This produces = 1.51(2), = 1.56(4), 
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FIG. 2. Temperature variations of resistivity and MR effect, (a) Resistivity measurements 
with current I along the a and c axes under applied magnetic fields of 0, 1 and 8 T. Ti = 100 
K labels the temperature where one structural anomaly occurs as shown in Fig. 1. The dashed 
lines are hts between 10 K and 40 K (details in text), and extrapolated to higher temperatures, 
(b) Corresponding MR values versus temperature. The MR effect along the a axis has a similar 
trend to that of the c axis albeit with a lower value. The positive MR values decrease sharply 
with increasing temperature below ~ 20 K, then gradually transfer into negative around Tn, and 
persist up to ~ 120 K. By contrast, in amorphous Gda;Sii_a; films, only the NMRV was observed 
[IMZ! near the metal-insulator transition. In (a), there is no big difference for the data at 0 and 
1 T above ~ 20 K. The solid lines in (a) and (b) in the dashed arrow direction in turn correspond 
to the ordinal axis-direction and applied magnetic-field strength as labeled, respectively. 


_ x. 43(2) and = 1.51(1). All w values are much smaller than 5 [6] indicative of 
anisotropic magnetic interactions. 

Anisotropic MR effect and positive MR value (PMRV) to NMRV transition. 

The most intriguing results from resistivity measurements are the anisotropic MR effect 
(Figs 2b, 3 and 4) and the existence of both positive and negative MR values, up to ~ 415% 
(comparable to the CMR value in manganites [HE] and one to two orders of magnitude larger 
than that of the RE-metals [2T]) and down to ~ -10.5% along the c axis at 8 T and 3 K and 
52.8 K, respectively. The MR anisotropy in the ac and be plans is shown in Fig. 3. They 
display a twofold symmetry at 7 K (Figs 3a and 3d). We notice that applied magnetic held of 
8 T does not suppress (produce) the (a) hump along the c and a axes, respectively, near Tn, 
and the MR twofold symmetry is persistent from 1 to 8 T, indicating that applied magnetic 
held in a strength of 8 T may just align or localize the 5d moments, and slightly rotate and tilt 
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FIG. 3. Angular-dependent MR values under different applied magnetic fields, (a) 

Angular-dependent MR values with applied magnetic field in a range of 0 to 8 T (1 T step) at 
7 K, and (b) at 120 K (0, 2 and 8 T). (c) For (a) and (b) measurements, current I is along the 
crystallographic c axis, and applied magnetic field, HoH, rotates away from the a axis with an angle 
of (p in the ac plane, (d) Angular-dependent MR values with applied magnetic field in a range of 
0 to 8 T (1 T step) at 7 K, and (e) at 120 K (0, 2 and 8 T). (f) For (d) and (e) measurements, 
/||a-axis and applied magnetic field, fJ-oH, rotates in the be plane with an angle of 7 deviated from 
the b axis. 


the 4/ moments while conserving the superzone energy gap. Therefore, the MR effect in GdSi 
exhibits a well separate feature of the temperature regions for the positive and the negative 
MR values (Fig. 2b), respectively, which is induced jointly by the AFM superzone effect [19] 
and the shift of the AFM transition to lower temperatures in external applied magnetic held 
analogous to the case in manganites muu- For metals, mean-held theories predict that 
spin huctuations induced by applied magnetic held from the localized magnetism produce a 
PMRV with the quadratic-held dependence in antiferromagnets, whereas a NMRV with the 
linear variation in ferromagnets and paramagnets [ 22 ] • However, in GdSi, the PMRV in the 
AFM state mainly displays a linear magnetic-held dependence (Figs 4a and 4c), while above 
Tui the absolute NMRV is proportional to the square of the strength of applied magnetic 
held (Figs 4b and 4d). Both the positive and negative MR ehects do not saturate at utilized 
maximum fioH = 8 T (Fig. 4). In addition, the ratio of the resistivity at 160 K and 7 K 
in Fig. 2a is already ~ 12-25, therefore, the cyclotron motion of the conduction electrons 
could be neglected at fioH = 8 T (i.e., UcT <C 1, where Uc is the cyclotron frequency and 
r is the life time of the conduction electrons). Therefore, these uncommon magnetic-held 
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FIG. 4. Field- and temperature- dependent MR values, (a) Field and temperature depen¬ 
dencies of the MR values with /||c-axis and applied magnetic field, along the a or c axis 

below Tji, and (b) above Tn. The representatives of the linear-field dependence of the PMRV (i.e., 
the PMRV is proportional to the strength of applied magnetic field) at 7 K (blow Tu) and the 
quadratic variation of the NMRV (i.e., the absolute NMRV is proportional to the square of the 
strength of applied magnetic field) at 60 K (above Tn) were shown in (a) and (b), respectively, (c) 
Field and temperature dependencies of the MR values with /||a-axis and applied magnetic field, 
fioH, along the 6 or c axis below Tu, and (d) above T^. 


variations indicate new transport mechanisms for the MR effects of GdSi. 

Discussion 

In GdSi, the condnction electrons (mainly 5d bands) are different from those responsible 
for the magnetism (4/ component pins possible part of the 5d component). The former is 
normally delocalized, acting as the magnetic glne among magnetic ions (Fig. 5), and scat¬ 
tered by them, leading to electrical resistance. The magnetism from the 4/ part is generally 
localized with weak interactions. Therefore, the LRO AFM state originates mainly from 
the isotropic RKKY interactions throngh condnction bands [23]. The interaction between 
localized moments, M^oc, and itinerant ones, rrtiti, can generate an extraordinarily large 
Zeeman splitting in the mean-held approximation [23| 

E = ■ H + 2J{H)r^±ti ■ (d^ioc), (1) 

where g* is the spectroscopic splitting factors for the carriers, /ib is the Bohr magneton, J 
is the effective exchange coefficient and (M^oc) is the averaged local moment in the regime 
of band electrons. Since r^±t± is nsnally small, and the second term conld be very large 
(e.g., amonnting to fractions of an eV in the LRO magnetic state of En-componnds [24]). 
in addition, J is strongly associated with applied magnetic held by virtne of modifying spin 
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FIG. 5. Schematic illustration of the spin states with and without applied magnetic 
field, fioH, in different temperature regimes, (a-c) At zero magnetic field, (d-f) At applied 
magnetic field of ^qH. When T > Ti, spin moments more or less rotate (f), depending on the 
strength of and the size of MA, from a pure PA state (that is strictly observing the Curie- 
Weiss law as shown in Supplementary Fig. S4) in (c). When Tjn <T < Ti, the short-range AFM 
spins that are attributed only to the local 4/ moments appear (b) accompanied by the generations 
of polarized itinerate 5d spins, based on the deviation of the unit-cell volume from the Griineisen 
model shown in Fig. lb, and possible small amount of localized 5d spins according to equation 
(1). Applied magnetic field mainly aligns the local AFM spins (e), leading to a decrease of the 
electron-spin scattering and resultant the NMRV. When T < T^ 2 , the LRO AFM state with almost 
equivalent AFM and FM interactions (see Supplementary Fig. S4) forms (a) with more itinerate Sd 
moments (based on the large PMV effect shown in Fig. lb). Applied magnetic field mainly localizes 
more f)d moments by enhancing the exchange of J in equation (1), resulting in the formation of 
magnetic polarons and the consequent PMRV (d). 


fluctuations of Mioc, the formation of magnetic polarons in the hd bands by this splitting in 
the LRO AFM state of GdSi is thus possible. Therefore, when T < Tn 2 , the modihed J at 
8 T drives some of the conducting moments (as foregoing remarks) to form local magnetic 
polarons that lead to a largely degenerate conduction (i.e., PMRV) (Figs 5a and 5b). In 
this case, the more extended 5d bands almost certainly offer a small FM component, which 
dominates the linear-magnetic-held dependence below Trj 2 . 

Above Tjii, the LRO AFM 4/ moments disappear, which is implied by the change of the 
X relationship between the a, h and c axes (see Supplementary Fig. S3a). However, the 
upward deviation of 1/x from the Curie-Weiss law below Ti (see Supplementary Fig. S4) 
along with the concomitant anomalous volume lattice distortions (Fig. 1) indicates that a 
kind of striking AFM state exists there, which can be ascribed only to the local short-range 
AFM 4/ moments (Fig. 5b) according to equation (1) (broad 5d states indeed could not 
form a magnetic state by themselves). This local AFM state theoretically [22] leads to the 
quadratic-magnetic-held variation (Figs 4b and 4d). The local AFM moments above Tm are 























possibly oriented irregularly in the PM state and become compulsorily aligned and to some 
extent spin-tilted (Fig. 5e), leading to the decrease in resistivity and thereby promoting a 
NMRV. 

Including the contribution of magnetic polarons into the mean-held approximation of the 
MR effect of AFM metals |22l l2SH27j, the MR value (MRV) then equals 

|MRV| 

«.4(|AlSP + 2|iiJS||iiJS|cosa+|iiJ“n, (2) 

where R is a constant, Min is the induced local magnetization by applied magnetic held, poM, 
a is the angle between Mj^ and M^^, and ^i^ is the magnetic susceptibility 

of the 4/ site along the AFM sublattice direction. At Tm < T < Ti, neglecting |Mj^|, 

I NMRV I oc |Mj^p = resulting in a quadric-magnetic-held dependence. At T < 

rK 2 , in our case, the linear part in equation (2) is dominate, i.e., PMRV oc |Mj^||M“| cos a = 

l^inll^inlho-ff cos a, which produces not only the linear-magnetic-held dependence but also 
the twofold MR symmetry (Figs 3a and 3d). Quantitative analysis requires the knowledge of 
the exact directions of ^i„ and M“. The involvement of magnetic polarons in understanding 
the MR ehects here is further supported by the existence of strong FM Gd-Gd interactions 
(see Supplementary Fig. S4). 

RE 4/ electrons generally remain localized so that their properties in an alloy closely 
resemble those in the single-elemental metals. It is therefore interesting to compare the 
microscopic conducting mechanism of GdSi with that of the 3d transition metals [2H1 1^ - 
Roughly, the 4/ shells highly screened by the 5s and 5p states are better shielded than 
the transition metal (TM) 3d ones. Therefore, the 4/ electrons are well embedded within 
the atom and the 5d and 6s states act as conduction electrons in the metals, though band 
structure calculations suggested that the upper part of the 4/ spin minority may hybridize 
with the unoccupied 5d and 6s bands On the contrary, for 3d TMs, the conductivity 
are mainly dominated by the s-s and d-d transitions at low temperatures, in addition to 
the s-d transitions (increasing the effective electron mass and shortening the mean free 
path) at high temperatures [281129] , supposing that both s and d electrons are conduction 
electrons. Since almost all RE 4/ metals and only some of the 3d TMs (Fe, Go, Ni and Mn) 
order magnetically at low temperatures, magnetism plays a more crucial role in controlling 
electrical property of RE 4/ metals. 

In summary, GdSi displays a wide array of novel behaviors including a two-step AFM 
transition, an anisotropic spontaneous MS effect, magnetostrictive strains, in particular, 
an exotic NTVE (spontaneous PMV) and a peculiar positive to negative MR transition, 
indicating a strong coupling between spin, lattice and charge degrees of freedom. The 
peculiar PMRV to NMRV transition and their nontrivial [22] magnetic-held dependencies are 
predominated by the nature of the 4/ moments and the polaronic 5d carriers. Goexistence 
of these extraordinary behaviors in the same Gd-compound is unique. The present results 
make GdSi an fascinating system for theoretical and further experimental explorations. 

Methods 

Since the neutral Gd is a strong thermal neutron absorber, we therefore performed a 
powder x-ray diffraction study of the pulverized GdSi single crystal on an in-house diffrac¬ 
tometer employing the copper Kai = 1.5406(9) A as the radiation with a 20 step size of 
0.005° from 10 to 300 K to explore the temperature-dependent structural modihcations. 
X-ray powder diffraction data were analyzed by Fullprof suite |T]. ZFG dc magnetization 
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(M) measurements were performed on a Quantum Design MPMS-7 SC quantum interfer¬ 
ence device (SQUID) magnetometer. The ZFC electrical resistivity (p) of bar-shaped (~ 
0.6 X 1 X 6 mm^) single crystals was measured by standard DC four-probe technique using 
a commercial physical property measurement system (PPMS), equipped with the option of 
rotating the sample in situ. Commercial silver paint and 25 pm gold wire were used for 
electrical contacts. MHV{poH,T,cj)/'y) = ^ X00%, where p{pqH,T, (p/'y) 

and p(0,T, 0 / 7 ) are the resistivity with and without pqH at a given temperature, T, and a 
rotating angle of 0 or 7 (Figs 3c and 3f). 
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Supplementary information for “Possible magnetic-polaron-switched positive and 
negative magnetoresistance in the GdSi single crystal” 

X-ray powder diffraction refinements and anisotropic stains. The observed and 
refined x-ray powder patterns at 10 K and 300 K are shown in Fig. Sla. The listed 
reliability factors in the caption indicate good refinements. In Fig. Sib, the full width at half 
maximum of the nuclear (1 1 2) Bragg peak by Lorentzian htting increases largely below T^, 
indicating the existence of strong strains that are produced by the gradual formation of the 
antiferromagnetic (AFM) state that is a process accompanied by the AFM sublattice domain 
movements. We first measured the standard LaBg powders, with which we extracted the 
instrumental resolution. Then we refined the collected x-ray powder patterns at 10 K with 
the obtained resolution using Fullprof suite [T] to extract the stain configurations as shown 
in Fig. S2. The existing AFM-driven anisotropic strains conserve the lattice symmetry 
or may lead to possible nuclear structural domain changes, which, however, may be too 
small to be detected by the instrument accuracy. It is stressed that the lattice change in 
GdSi, e.g., = 3.2(1) x 10““^, is much smaller than that of the symmetry-broken 

AFM metals such as SrFe 2 As 2 , where = 3.1(2) x 10“^ [2]. Similar AFM-driven 

strains were also observed in the classical 3d Chromium metal that also displays charge- 
and AFM-spin-density waves [3]. 

Magnetization measnrements, AFM transitions and magnetic anisotropy. To 

conhrm that the positive to negative magnetoresistance (MR) transition is really induced by 
the shift of Tn to low temperatures with applied magnetic helds, we exactly hxed the crys¬ 
tallographic a and c axes, taking advantage of the rotator, and measured the temperature- 
dependent magnetization as shown in Fig. S3. It is unexpected that the AFM transition 
experiences two steps at Tm and Tu 2 (Fig. S3a), respectively. From Tm to Tk2, there is 
a slight decrease in magnetization, whereas below Tu 2 it hrst-order-like diminishes rapidly 
indicative of the formation of the long-range-ordered AFM state. Note that Xb > Xc > Xa 
below T[j2 transfers to Xc > Xb > Xa between Trj2 and Tm, and subsequent to Xc > Xa > Xb 
above Tm, indicating the existence of magnetic anisotropy (MA) and different magnetic 
states in the three temperature regimes. The observed MA here is in agreement with the 
anomalous lattice constants (Fig. la), the lower w values than 5 (Fig. 2a), and the ap¬ 
proximately zero paramagnetic Curie temperatures as discussed below (Fig. S4). Applied 
magnetic fields indeed shift Tm and Tu 2 to lower temperatures in the constructed HqH-T 
phase diagram shown in Fig. S3b indicative of some local Cd-Cd AFM interactions, which 
is dramatically different with the purely itinerant Cr whose is strongly rigid in response 
to applied magnetic fields even up to 14 T, nevertheless extremely sensitive to strains [3]. 

We also measured the temperature dependence of the magnetic susceptibility y = M/hqH 
along the a and c axes at 0.06 T in a temperature range of 3-300 K. The corresponding 
inverse susceptibility 1/y is shown in Fig. S4. Fitting the data between 150 K and 300 K 
to the Curie-Weiss (CW) law, we obtained the effective magnetic moment, /Xeff, and the 
paramagnetic Curie temperature, 6*p, to be = 8.437(2) /iB/Cd and = 0.33(1) K 
for hqH along the [100] direction, and = 8.424(1) /is/Cd and = 0.09(7) K for hqH 
parallel to the [001] direction. Below ~ 100 K, y“^(T) upward deviates largely from the 
CW law, implying an increasing contribution of the AFM spin correlations. It is pointed 
out that small upward deviation already appears at ~ 120 K consistent with the observed 
negative MR effect at 120 K as shown in Figs 3b, 3e, 4b and 4d. However, this kind of small 
local AFM spin states do not generate appreciable anomalies in the structural parameters 
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FIG. 6. Powder x-ray diffraction patterns and refinements, a, Observed (circles) and 
calculated (solid lines) x-ray powder-diffraction patterns at 10 K (down) and 300 K (up). The 
vertical bars mark the positions of Bragg reflections. The lower curve represents the dif¬ 
ference between the observed and calculated patterns. The reliability factors of refinements 
are: 10 K [Rp{%) = 1.06,i?^p(%) = 1 . 57 , = 1.2,7, Rb{7o) = 2.65,i?F(%) = 2.84]; 300 K 

[Rp{7o) = 0.961, i2^p(%) = 1.33, = 1.22,Rb{7o) = 1.47, i2i7’(%) = 1.70]. b, Observed nuclear 

(1 1 2) Bragg peaks at selected temperatures from 10 K to 70 K (below T ^2 and above Tjn). The 
solid lines are fits of the Lorentzian line-shape, and the dashed line is a guide to the eye. The small 
horizontal bar indicates the instrumental resolution. The fitted peak position and full width at 
half maximum are labeled at respective temperatures. 



FIG. 7. Strain patterns of GdSi. Extracted anisotropic stain configurations in the ac (a), ba 
(b) and be (c) planes at 10 K. 


until upon cooling to a temperature of ~ 100 K, where a sufficiently large magnetoelastic 
coupling that is able to shake the crystalline lattice may be just attained, as shown in Fig. 1. 
The occurrence of very small positive paramagnetic curie temperatures, 6p, in the dominate 
AFM matrix, indicates that either the 4/ moments in GdSi are purely localized without any 
interaction at all, or there exist mixed ferromagnetic (FM) and AFM Gd-Gd interactions 
with an approximatively counterbalance of their strengths like the rare case in AFM ZnGr 2 S 4 
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FIG. 8. Antiferromagnetic transitions and ^qH-T phase diagram, a, Temperature- 
dependent magnetization under ^jlqH = 0.06 T (along the a, b and c axes) and 2.5 T (along 
the a and c axes, shifted down for comparison). The solid lines are linear fits of the a-axis mag¬ 
netization (0.06 T) between Tii 2 (~ 52.94 K) and Tm (~ 54.3 K), and just below Tn 2 (50-52.94 K), 
indicating a clear slop change from 7.4(1) x 10“® (above Tji 2 ) to 8.8(1) x 10“® /tb/K (below Tji 2 ). 
b, Measured fj-oH-T phase diagram around Tm and Tn 2 . Error bars are the standard deviation 
based on entire measurements. The solids are guides to the eye. The intermediate (I) phase bor¬ 
dering the long-range-ordered static 4/ AFM state and the local short-range 4/ AFM spins (AS) 
may be attributed to a short-range-ordered 4/ AFM state because of the large suppression of Tm 
and Tji 2 , consistent with the stiffness of the hump in resistivity under applied magnetic fields (Fig. 
2a). 




FIG. 9. Inverse susceptibility versus temperature. As measured (circles) at 0.06 T along the 
a (a) and c (b) axes in GdSi. The solid lines indicate a Curie-Weiss behavior between 150 K and 
300 K, which was extrapolated to zero temperature (dashed lines). One anomaly displays at ~ 5 
K. Ti = 100 K labels the temperature where one structural anomaly occurs as shown in Fig. la. 

[1] . The coexistence of FM and AFM exchange interactions often causes spin frustrations. 

Magnetic anomaly at ~ 5 K. The anomaly at ~ 5 K in Fig. S4 may suggest a different 
magnetic state below 5 K with the one between 5 K and Tji 2 . A similar anomaly was also 
observed in Ref. [5]. The htting of the resistivity between 10 K and 40 K (Fig. 2a) also 
produces = -7.3(2) x 10“^ mO mm, = 1.5(1) x 10“^ mO mm, = -1.0(2) 
X 10“^ mO mm and = 1-1(1) X 10“^ mO mm. It is interesting that both pg°^ and 
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Pg’°^ (at zero K) are negative, indicating that at lower temperatures the conduction mode 
may be modihed j6], or there exist magnetic transitions consistent with the observed clear 
anomalies at ~ 5 K in magnetization measurements (Fig. S4). 

We did not observe any anomaly indicative of the weak FM order at ~ 78 K shown in 

Ref. [7]. 

Future studies. Determining the magnetic structures and interactions will be the sub¬ 
jects for future studies, which may be not easy because of the strong thermal neutron 
absorption of neutral Gd and the intangible influence of 5d conducting moments on the 4/ 
AFM ones. 
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